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Interfacial phase change memory (iPCM), that has a structure of a superlattice made of alternating 
atomically thin GeTe and Sb 2 Te 3 layers, has recently attracted attention not only due to its superior 
performance compared to the alloy of the same average composition in terms of energy consumption but 
also due to its strong response to an external magnetic field (giant magnetoresistance) that has been 
speculated to arise from switching between topological insulator (RESET) and normal insulator (SET) 
phases. Here we report magneto-optical Kerr rotation loops in the visible range, that have mirror symmetric 
resonances with respect to the magnetic field polarity at temperatures above 380 K when the material is in 
the SET phase that has Kramers-pairs in spin-split bands. We further found that this threshold temperature 
may be controlled if the sample was cooled in a magnetic field. The observed results open new possibilities 
for use of iPCM beyond phase-change memory applications. 



Topological insulators (TIs) have attracted increasing attention from both solid-state physics and applications 
in spintronics devices perspectives 1 " 5 . TIs such as Sb 2 Te 3 and related chalcogenides built from quintuple 
layers (QL) satisfy both spatial inversion symmetry and time reversal symmetry, resulting in a single Dirac 
cone at the T point 6 ' 7 , where the spin bands are degenerate due to the presence of inversion symmetry. However, 
once spatial symmetry is broken, e.g. by an external electric field, the band gap opens and the spin degeneracy is 
lifted resulting in a Rashba-split phase 8 . To date, experimental studies of TI properties have usually been carried 
out using single crystal surfaces in vacuum 4 ' 6 . 

Recently, it was predicted that interfacial phase change memory (iPCM) 9 , which consists of a superlattice with 
alternating Sb 2 Te 3 (TI) and normal insulating (NL) GeTe spacer layers: [(GeTe) 2 (Sb 2 Te 3 ) i ] n , where / and n are 
integers 10 , originally developed with the objective to reduce energy consumption, may exist in different topo- 
logical phases 11 . In particular, it was argued that iPCM with / = 1 was a Dirac semimetal with bulk gapless states 12 . 
Furthermore, electric-field induced giant magnetoresistance was experimentally observed at room temperature 
despite no magnetic elements being included 13 . It has been speculated that this phenomena is caused by the 
Rashba effect due to the breaking of spatial inversion symmetry 12 . 

In this work we report on optical detection of topological properties of iPCM, in particular, on the use of the 
magneto -optical Kerr rotation. Magneto -optical phenomena such as Kerr- or Faraday- rotation were proposed to 
be alternative tools to probe the spin-dependent response of TIs, besides the spin transport analysis 1415 . The rather 
simple polar-Kerr-rotation technique, where linear polarized light and an external magnetic field are both applied 
normal to a sample surface, is especially appealing 16 . The light induces an AC electrical field in the illuminated 
spot along the polarization direction, generating a time dependent Hall conductance t), which affects the off- 
diagonal elements of the dielectric tensor (i, j = x, y, z) as the dielectric constant arises from the electrical 
conductance tr^ within the framework of Kubo's linear response theory 17 . Recently, Faraday and Kerr-rotation on 
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TI surfaces were theoretically discussed 18 , and it was predicted that a 
resonant condition may induce a giant Kerr rotation d K by — nil 
radian in the low frequency (terahertz) range for a wavelength far 
longer than the TI film thickness 18 . Here we demonstrate that this 
approach is not limited to the resonance between the broken Dirac 
cone bands that form the band gap between the filled and empty 
states but even using visible light, where absorption involves transi- 
tions between the valence spin bands and higher conduction bands, a 
large Kerr-rotation can be observed. 

Results 

SET&RESET structures of interfacial phase change memory 
(iPCM). Figure 1(a) shows two iPCM phases (the high-conductivity 
SET and low-conductivity RESET phases), the corresponding 
electronic band structures calculated including spin-orbit coupling 
(SOC) effects in bulk 12 (See Method), and the high-angle annular 
dark field transmission electron microscope (HAADF-TEM) images 
of the RESET and SET phases of the iPCM [(GeTe) 2 (Sb 2 Te 3 ) 1 ] 20 
structure fabricated at 520 K on a Si substrate. The RESET phase 
has spatial inversion symmetry with a Dirac cone in the bulk while 
the SET phase doesn't. The SET and RESET phases reversibly 
transform into each other within a temperature range of 380 K and 
450 K as shown in Figure 1(b), where the temperature dependence of 
the resistance of iPCM is shown. The shown behavior is highly 
unusual for phase-change memory where the change of resistance 
upon cooling the SET phase is typically negligible 19 . It should also 
be noted, that iPCM in the SET state has an unusually structured 



optical absorption spectrum in the visible wavelength range as 
shown in Fig. 1(c). In particular, despite a rather small band gap of 
—0.3 eV, the SET phase exhibits a sharp absorption increase near 
700-nm, which corresponds to transitions from below the Fermi 
level to unoccupied energy bands located ca. 1.8 eV (690 nm) above 
the Fermi level (see the band structure of the SET phase in Fig. 1(a)). 
In contrast with the band structure in the RESET phase, the spin 
degeneracy is lifted in the SET phase by the breaking of spatial 
inversion symmetry generating a narrow band gap, while the phase 
still preserves the time reversal symmetry. 

Mirror-symmetric Kerr-rotation of iPCM thin films. Polar Kerr- 
rotation results of the RESET and SET phases of [(GeTe)2(Sb 2 Te 3 ) 1 ]2o 
iPCM obtained using 690-nm light are shown in Fig. 2 (a), measured 
at room temperature and at 440 K (See Method). While the signal at 
room temperature was within noise, a strong and unusual Kerr-loop 
with 0^ ~ —0.3 degree was observed at 440 K within the applied 
magnetic field range (±13 kOe). The loop has four resonances with a 
mirror- symmetry completely dissimilar to the deformed-rectangular 
hysteresis loop with respect to the magnetic field polarity typically 
observed in ferromagnetic materials. The drawing order of the shown 
curve is as follows: 0 ^ 1 ^ 2 ^ 3 ^ 4 ^ 5(0) ^ 6 ^ 7 ^ 8 ^ 9 ^ 
0. The magnetic field required to generate resonance peaks 1 and 3 
monotonically increased with decreasing temperature as shown in 
Fig. 2(b). Once the temperature was decreased to less than the 
threshold (—380 K), the signals disappeared abruptly, but re- 
appeared reversibly when the temperature was increased again 




(a) (c) 

Figure 1 | (a) Two different crystalline phases (RESET and SET) ofiPCM [(GeTe)2(Sb 2 Te 3 ) 1 ] n , the corresponding band structures (only shown in M-r-K 
k-space) calculated by ab-initio computer simulations, and high-angle annular dark field transmission electron microscope (HAADF-TEM) images of the 
real RESET and SET phases with superimposed simulation models. Sb 2 Te 3 -quintuple layers and (GeTe) 2 layers were clearly identified as well as a vacancy 
layer, respectively. The spin-bands are degenerate in the RESET phase while the degeneracy is lifted in the SET phase, (b) Temperature dependence of the 
resistance upon heating (red) & cooling (blue) demonstrating a reversible change, using a [(GeTe) 2 (Sb 2 Te 3 ) 1 ] n = 20 superlattice film, (c) Optical 
absorption of the SET phase calculated using ab-initio computer simulations (blue), and the zoom of the visible wavelength range (inset) with an 
experimentally observed curve using a [(GeTe) 2 (Sb 2 Te 3 )!] n = 2 superlattice film. 
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Figure 2 | (a) Mirror- symmetric Kerr-loop observed from iPCM [(GeTe)2(Sb 2 Te 3 ) 1 ]2o structure at 440 K (blue). For comparison the result of the room 
temperature measurement is shown in black. The order of the Kerr signal with variation in magnetic field was: 0 — » 1^2^3^-4— >5(0)^6^7^-8 

9 ^> 0; (b) Magnetic fields corresponding to 1 and 3 resonances as a function of temperature; below 380 K the Kerr rotation disappears; (c) The 
Kerr-loop arising from an 1-nm thick (GeTe) 2 layer on top of a 5-nm Sb 2 Te 3 layer (used to orient GeTe) at 440 K; (d) Kerr-signal behavior of (c) due to a 
magnetic field at 440 K (red) and by cooling the sample temperature at 1 1.85 kOe (blue). 0 to 9 and arrows marked on the loop in (a) are the order of one 
Kerr-loop cycle. 



above the threshold. The threshold temperature for the (dis) appea- 
rance of the Kerr-signal is in good agreement with the transition 
temperature between the SET and RESET phases observed in the 
electrical resistance measurements (Fig. 1(b)), suggesting that the 
observed Kerr loops are associated with the SET phase. 

The sense of the Kerr rotation changed randomly for various 
samples and even different pieces of nominally the same sample 
but the mirror symmetry and the threshold nature of the Kerr rota- 
tion were reproduced for all samples. We further found that the 
signal intensity dropped significantly when the light wavelength 
was changed from 690-nm to 650-nm (Fig. SI). At the same time 
the magnetic field at which the resonant Kerr rotation occurred 
shifted to a lower value for the light with a higher photon energy 
suggesting that the photon energy and Zeeman splitting add up to 
generate the necessary conditions for the Kerr rotation resonance. 
The Kerr rotation angle became negligibly small for 400-nm light, 
further demonstrating that the process has a resonant nature. 

To study the origin of the magneto -optical response in iPCM, we 
also examined the phenomena using several different films. As 
shown in Fig. 2(c), a 1-nm thick GeTe layer fabricated on a 5-nm 



thick Sb 2 Te 3 layer (required to orient the growth along the < 1 1 1 > 
direction) showed similar loops at 440 K. At the same time, the Sb 2 Te 3 
layer alone, and a 20-nm thick epitaxial GeTe layer grown on Si (1 1 1) 
did not exhibit any measurable signal. This result underscores the 
important role of GeTe/Sb 2 Te 3 interfaces in the observed effect, sim- 
ilar to the case of the giant magnetoresistance observed in iPCM but 
not in alloyed materials of the same average composition 13 . It is worth 
noting that the two resonance drops for the GeTe layer shown in 
Fig. 2(c) merge into one demonstrating a rather small hysteresis, 
compared to the much larger hysteresis in [(GeTe) 2 (Sb 2 Te 3 ) 1 ] 20 , a 
result that suggests that the superlattice is ferro-magnetized and has 
the potential to store spin in the heterostructure due to proximity- 
induced interfacial magnetization 20 ' 21 . 

It is interesting to note that, when the 1-nm thick GeTe layer on the 
5-nm thick Sb 2 Te 3 layer was cooled under the applied magnetic field, 
the Kerr loop could be observed down to 320 K, i.e. 60 degrees lower 
that the threshold temperature without the magnetic field (Fig. 2(d)); 
the same phenomenon was observed in [(GeTe) 2 (Sb 2 Te 3 ) 1 ] 20 down to 
310 K as well (not shown), i.e. the phase transition between SET and 
RESET could be controlled by the external magnetic field. 
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Magneto-electric effect in the RESET phase. Although the RESET 
phase does not exhibit Kerr rotation, it also has a specific response to 
the applied magnetic field. As shown in Fig. 3, the optical dielectric 
constant of the [(GeTe)2(Sb 2 Te 3 ) 1 ]2o film at room temperature was 
increased in the region, where the magnetic field (4 kOe) was 
applied. In the Maxwell equations, a static magnetic field does not 
couple to a dielectric constant. This phenomenon is highly unusual, 
except for a few multiferroic composites with highly-oriented 
interfacial contact of ferromagnetic and ferroelectric thin films 22,23 . 

Discussions 

We now address a possible mechanism for the observed phenom- 
enon. While the accurate description of the observed results requires 
further study, the basic results can be understood as follows. As 
shown in Fig. 1(a), the degeneracy of the spin-bands in the SET phase 
is lifted due to the breaking of spatial inversion symmetry by the 
phase transition from RESET to SET states. 

Before we discuss the origin of the Kerr-loop in the SET phase, we 
consider the system in the RESET phase that possesses both time 
reversal and spatial inversion symmetries. In the absence of an 
external magnetic field, the up -spin ( j) and down -spin ( j) electrons 
satisfy the condition E(—k, [) = E(+k, |), where E and k are band 
energy and electron momentum, respectively. Since in a two-dimen- 
sional system with the Dirac cone spins are usually oriented in 
plane 24 , it may be more appropriate to re-write the above condition 
as E(—k, = E(+k, -^). Because the energy of the spin-orbit 
interaction is oc {(7 x k y — (7 y k x ) E z , where g x and a y are Pauli spin 
matrices, and E z is the z-component of the electrical field induced by 
the laser beam (the z-direction is defined to be normal to the sample 
surface), it is essentially impossible for our laser set-up with just four 
degree away from normal incidence to induce a strong E z lifting of 
the band degeneracy in the RESET phase shown in Fig. 1(a). Even for 
the case that the actual film surface is not perfectly normal to E& 
Kerr- rotation components 0^+ and 0^- cancel each other due to the 
presence of degeneracy. 

When time reversal symmetry is broken by an application of an 
external magnetic field B z (] ) ex or B z ([ ) ex , the situation becomes dif- 
ferent due to the Zeeman effect. The Zeeman effect in a thin film 
topological insulator was recently theoretically investigated for the 
case when an external magnetic field was applied normal to the 
surface 25 ' 26 . The magnitude depends on the hybridization energy of 
the helical spins between the top and bottom surfaces of the thin film 
topological insulator. Although the Zeeman energy Az = 1/2- g/n B B is 
usually small, it is noted here that in the presence of large spin-orbit 
coupling, the Lande factor g of topological insulators can be rather 



large, over 50 for for Bi 2 Se 3 , which corresponds to Az of a few tens 
me y27,28 gj nce the helical spin coupling between the top and bottom 
interfaces in the superlattice and spin-orbit coupling give rise to a 
Dirac cone in the bulk state in the [(GeTe) 2 (Sb 2 Te 3 ) 1 ] 2 o iPCM film 12 , 
a large Zeeman effect can also be expected by analogy in our case. The 
large change between the band structures originates from the helical 
spin coupling. In addition, it was theoretically reported that a super- 
lattice structure built from [(BiTeI)/(Bi 2 Te 3 ) m ] blocks, where / and m 
are integers, also has a large Rashba splitting due to the breaking of 
spatial inversion symmetry, which leads to an out-of-plane spin 
component (z) 29 . The Zeeman energy shift may thus be comparable 
to the band split (—50 meV) near 1.8 eV corresponding to the laser 
wavelength (690 nm) as shown in the SET phase (Fig. 1(a)). The 
magneto -optical effect shown in Fig. 3 is presumably attributed to 
the Zeeman effect in the iPCM film. 

Due to the combination of the Zeeman effect and the breaking of 
time reversal symmetry, two different cases can be considered 
depending on the direction of the external magnetic fields: B z (] ) ex 
and B z (l) ex . As the internal magnetic field from the spin-orbit coup- 
ling is added to the B z (\or[ ) ex in the superlattice, hereafter we use the 
effective magnetic field, B z (]ox[) e ff. For increasing B z (]) ex , the 
energy balance between spin-up and spin-down bands is changed 
to E(-k, l)>E(+k,] ), as shown in Fig. 4. (a) and (b). While B z (] ) eff 
is within a narrow resonant region, the balance between the clock- 
wise d K+ and anti- clockwise 0^_ Kerr-rotation components is bro- 
ken, with the former dominating the process. It is notable that the 
energy difference between the spin-up and spin-down sub-bands is 
small as shown in Fig. 1(a) -SET, and is mostly induced by the loss of 
spatial inversion symmetry. Therefore, the resonance condition is 
limited to a narrow B z (]) e ff range (Fig. 4(b)). As a result, a steep 
Kerr-rotation appears. When B z (]) e jj is further increased beyond 
the resonant region, the efficiency of both d K+ and 0^_ transitions 
decreases, resulting in disappearance of the signal (Fig. 4 (c)). 

When the magnitude of the magnetic field is decreased, the system 
undergoes the reverse change. However, once the SET phase is mag- 
netized by a combination of the external magnetic field, the return 
path (2^5) is different 30 : the magnetic moment M z (\) due to the 
magnetization must be added to B z (] ) e ff. Therefore, when B z (] ) ex is 
decreased from "2" to "3", the resonance condition is not attained at 
the same B z (]) ex as for the outward path (position 1), resulting in 
hysteresis at position 3 (Fig. 4 (d)). From 3 to 5, the residual mag- 
netization is gradually lost through scattering processes when B z (\ ) is 
decreased. Finally, the spin (up and down) density balance recovers 
alongside with time reversal symmetry at B z = 0. As a result, along 
the return path the Kerr-rotation angle has a peak at lower applied 
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Figure 3 | Dielectric constant obtained from ellipsometric measurements using a wavelength of 633-nm. A 10 X 10 mm 2 permanent magnet (4 kOe) 
was placed under the sample (left), and the dielectric constant at 633-nm was measured across the sample including areas above the magnet and on both 
sides (right). Blue and red curves represent the magnetic polarities N and S pointing up, respectively. The color section in light blue is the region, which the 
magnet was placed. 
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Figure 4 | The basic Kerr-loop mechanism of the SET phase, (a) Splitting into two spin-bands in the presence of time reversal symmetry in the SET phase 
at zero magnetic field B zex = 0. The clockwise and anti-clockwise Kerr-rotation components cancel each other, (b) Resonance state due to E(+k x , j) at 
B z (\) resonance — B z (\) e ff. (c) Non-resonance state when B z (\) e ff> B z (]) resonance . (d) Resonance state on the return path. As the SET phase is once 
magnetized with M z (]) on the outward path, the return resonance state appears at a lower external magnetic field than B z (]) ex on the outward path. The 
details depend on the band structure and the value of the magnetic field applied, as illustrated by the four panels. The notations 0,1,2,3, and 4 correspond 
to the positions marked with the same indices in Fig 2(a). 



magnetic field and decreases less abruptly compared with the onward 
path (Fig. 4 (a)). 

We now turn to the important issue of the mirror symmetry of the 
Kerr-loop. The off-diagonal elements in the dielectric tensor are 
proportional to the transition probabilities of the electrical dipole 
polarization operators P+ or P_ between energy bands, and the 
Kerr- and Faraday- rotation can be calculated from the summation 
A of individual oscillator strengths, (f±) mn = oc-oo mn |P ± | 2 , i.e. 
A=2?{(f+) mn — (f-) mn }> where m, n, oo mn and a are the band labels, 
the resonance frequency between m-th and n-th bands, and a con- 
stant, respectively 31 . Generally in magnetic materials, the sign of A 
does not reverse with the polarity of the external magnetic field. As a 
result, a normal magnetic hysteresis loop is generated from the off- 
diagonal elements, s } 



yx 



— 8 



xy 



We note that if the Kerr rotation in the SET phase is driven by 
broken time-reversal symmetry due to the magnetic field, the Kerr 
rotation should change sense when the magnetic field is flipped. It 
suggests that the present Kerr signal is due to the presence of gyro- 
tropic order in the iPCM-SET phase, which has recently been dis- 
cussed for cuprate superconductors 32,33 and lT-TiSe 2 34 ; in the latter it 
was associated with charge density waves. Gyro tropic order breaks 
mirror symmetry, which causes optical activity. In the present case, 
mirror symmetries with respect to any plane perpendicular to the 
surface should be absent, because such symmetry makes the Kerr 
rotation signal for normal incidence to be an odd function of the 
magnetic field. 

The next question is where the gyrotropic order originates from. 
The facts that the sense of the Kerr rotation in the SET phase 
randomly varies among different nominally similar samples, while 



for each sample the signal is perfectly reproduced when the sample is 
cooled down to the RESET phase and warmed up again to the SET 
phase suggests that the gyrotropic order (handedness) is determined 
by the superlattice structure itself, rather than by spontaneous order- 
ing of electronic states. We associate the presence of gyrotropic order 
with the sequence of Ge and Te planes in the GeTe block of the SET 
phase. Indeed, while the rhombohedral structure of GeTe is often 
described as distorted cubic, it should be noted that there is a pro- 
nounced bond length (energy) hierarchy with the shorter Ge-Te 
bonds of 2.80 A and longer Ge-Te bonds of 3.13 A. Consequently, 
in a hexagonal representation the structure can be seen as layered 
with covalently bonded layers held together by weaker van-der- 
Waals interaction (Fig. 5(a)). The unit cell contains three covalently 
bonded layers with the Ge (marked using different colors for differ- 
ent layers) positions varying among the three layers. Shown in 
Fig. 5(b) is a GeTe unit cell with additional Ge atoms along the c 
axis. In the figure, the Ge atoms are connected by "bonds" that serve 
to guide the eye. One can clearly see that the Ge atoms form a helix 
oriented along the c-axis, i.e. along the propagation direction of the 
light. The Ge helix is further visualized in Fig. 5(c), where a view from 
a different angle is presented. Based on this model, the mirror sym- 
metric Kerr rotation is expected in the SET phase. 

Spatial inversion symmetry of the RESET phase (and of its GeTe 
block in particular) prevents the formation of helices in the latter. It 
should also be noted that while the existence of the mirror symmetric 
Kerr rotation per se arises from the presence of helical fragments, its 
resonant nature is a result of the Zeeman splitting by the magnetic 
field and the breaking of the spatial inversion symmetry; with both 
the helical fragments and the resonant band splitting being essential 
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Figure 5 | (a) GeTe in a hexagonal representation. Te atoms are shown in 
tan and three inequivalent Ge atoms are shown in red, green, and blue. The 
shorter Ge-Te bonds within covalently bonded layers are shown as dual- 
band cylinders and the unit cell is shown as black lines, (b) The unit cell of 
GeTe with several additional Ge atoms along the c-axis. The Ge atoms are 
connected by gray lines to guide the eye. One can see that the Ge atoms 
form a helix, (c) A view on the Ge helix from a different angle. Only the Ge 
atoms and the unit cell are shown. Te atoms are invisible to make the 
helical structure unobscured. 

attributes of the SET phase. Further work is underway to elucidate 
this issue in more details. 

In conclusion, although iPCM superlattices were first designed for 
reduce the switching energy in memory devices, the present work 
demonstrates that they may become an interesting platform for 
experimental investigations of topological properties in engineered 
structures. The reversible band gap control by an external electric 
field in the RESET phase, and the spin- split band in the SET phase 
existing in iPCMs at different temperatures are highly attractive for 
applications beyond nonvolatile memory. 

Methods 

Fabrication of iPCM films. iPCM films were fabricated on Si wafers (13—17 Q-cm) 
using a helicon- wave sputtering system that features a large (200 mm) target- 
substrate separation. The individual layers composing the (GeTe) x /(Sb 2 Te 3 ) y 
structure were fabricated using GeTe and Sb 2 Te 3 composite targets (2-inch) using an 
automated shutter control system with a substrate temperature of more than 500 K 
and pressures less than 0.5 Pa Ar. The unit thickness of the (GeTe) and (Sb 2 Te 3 ) 
layers were 0.4 nm and 1.0 nm, respectively. To ensure strong crystalline orientation, 
a 5-nm thick Sb 2 Te 3 layer was deposited prior to the fabrication of iPCM films. As a 
consequence, all the sub-layers were crystalline in the as-deposited state and exhibited 
a strong preferred crystallographic orientation towards the < 1 1 1 > direction. Finally, 
a 20-nm thick ZnS-Si0 2 layer was deposited as a cap layer to protect the films from 
oxidation. 



Measurements. The film resistance was measured in plane using a two-probe system 
equipped with a small heater-stage. As electrodes, 50-nm thick W films were 
deposited at the sample edges. Most of the Kerr rotation measurements were 
performed using a linear polarized laser beam, with 690- nm wavelength & 500 um 
spot size, incident at an angle of 4 degrees from the normal to the sample surface. The 
sample was set on a heater stage surrounded by a magnetic coil (maximum field: 
±13 kOe). The initial magnetic field was zero and the magnetic field change was in 
the following sequence: 0 kOe — > + 13 kOe — > 0 kOe^> -13 kOe^> 0 kOe. Since this 
system did not allow us to change the wavelength, we used a different system for 
measurements using 650-nm and 400-nm light. In the latter the magnetic field was 
swept not from zero but from the maximum positive to maximum negative values and 
back, namely, + 10 kOe -> 0 kOe -» - 10 kOe -> 0 kOe -» + 10 kOe. The optical 
dielectric constant was measured by an ellipsometer equipped with a He-Ne laser 
(633-nm wavelength). The magnetic field was applied to the sample using a 1 mm 
thick flat square permanent magnet (10 mm X 10 mm area) with a 4 kOe field. The 
laser beam was scanned every 2 mm to measure the optical parameters in between the 
sample edges. 

Simulation of the SET&RESET models and the corresponding band structures. 

The electronic band structures of the iPCM films were simulated using two ab-initio 
simulation codes: CASTEP and WIEN2K. IPCM models were first relaxed using the 
plane wave code CASTEP with a GGA exchange correlation term as constructed by 
Perdew-Burke-Ernzerhof (PBE). A 4 X 4 X 1 Monkhorst-Pack grid was used for 
integration and ultrasoft pseudopotentials were used with a cutoff energy of 
230 eV 35 ' 36 . After 0 K Broyden geometrical optimization, the models were transferred 
to WIEN2K 37 to include spin-orbit coupling effects. Wien2K is an all electron code 
that uses a linearized augmented plane wave + local orbital (LAPW + lo) basis within 
density- functional theory. The PBE exchange -correlation was used. A Monkhorst- 
Pack grid of 4 X 4 X 1 was used for integrations in the Brillouin zone and a R x K max = 
7.0 value was used for the plane wave component of the plane-wave basis used 
between augmentation spheres. 
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